Analysis of Boundary Slip in a Flow with an Oscillating Wall 
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Molecular dynamic (MD) simulation is used to study slip at the fluid-solid boundary in an un- 
steady flow based on the Stokes second problem. An increase in slip is observed in comparison to 
the steady flow for shear rates below the critical shear rate of the corresponding steady flow. This 
increased slip is attributed to fluid inertial forces not represented in a steady flow. An unsteady 
mathematical model for slip is established, which estimates the increment in slip at the boundary. 
The model shows that slip is also dependent on acceleration in addition to the shear rate of fluid at 
the wall. By writing acceleration in terms of shear rate, it is shown that slip at the wall in unsteady 
flows is governed by the gradient of shear rate and shear rate of the fluid. Non-dimensionalizing the 
model gives a universal curve which can be used to find the slip boundary condition at the fluid-solid 
interface based on the information of shear rate and gradient of shear rate of the fluid. A governing 
non-dimensional number, deflned as the ratio of phase speed to speed of sound, is identified to help 
in explaining the mechanism responsible for the transition of slip boundary condition from finite 
to a perfect slip and determining when this would occur. Phase lag in fluid velocity relative to 
wall is observed. The lag increases with decreasing time period of wall oscillation and increasing 
wall hydrophobicity. The phenomenon of hysteresis is seen when looking into the variation of slip 
velocity as a function of wall velocity and slip velocity as a function of fluid shear rate. The cause 
for hysteresis is attributed to the unsteady inertial forces of the fluid. 

PACS numbers: 47.11. Mn,47.61.-k,68.08.-p,83.50.Rp 



I. INTRODUCTION 

The no-slip boundary condition at the interface of a 
fluid and a solid wall has been the subject of many inves- 
tigation for more than a century [lH3 ■ Navier Q was the 
first to introduce the linear boundary condition, which 
was later also proposed by Maxwell [5|, and it remains 
the standard characterization of slip even today. Slip at 
the boundary although prevalent, is negligibly small in 
most continuum and macro-scale applications. Hence, 
the no-slip boundary condition is widely accepted and 
has been shown to give accurate results in such applica- 
tions. However, in many micro- and/or nano-scale ap- 
plications the first breakdown of continuum assumption 
often occurs at a solid boundary in the form of velocity 
slip. 

As transport in an ever smaller scales are considered, 
surface forces and effects start playing a more profound 
role on fiuid transport than the bulk forces. Results 
from various computer simulations (6l-[l7j , which has been 
backed up by a nunrber of laboratory experiments [T8l - [25j 
show the presence of slip in fluids at the boundary. The 
advent of molecular dynamic simulations proved to be 
a considerable aid in understanding slip, as performing 
experiments at such scales is difficult. Most MD simula- 
tions have been focused on steady flows. MD simulations 
of a shear-driven steady flow by Thompson and Robbins 
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|8[ showed dependence of slip on wall-fluid properties 
such as density of the wall relative to fluid, the strength 
of liquid-solid coupling, and the thermal roughness of the 
interface. They also observed layering of fluid normal to 
solid walls Q . Thompson and Troian performed MD 
simulation of the steady Couette flow. They observed 
that at small shear rates the boundary condition is con- 
sistent with the Navier model. However, as the shear 
rate is increased the Navier condition breaks down and 
the slip length increases rapidly with shear rate. They 
discovered that as the wall velocity is increased the slip 
at the wall is nonlinearly increased to infinity at a critical 
shear rate at which point the wall is no longer able to im- 
part any further momentum to the fiuid. They also went 
on to find a universal curve that gives the slip length for a 
specific shear rate irrespective of wal- fluid properties 0] . 

Boundary slip has been the subject of less investiga- 
tion in unsteady flows. A few exceptions are in unsteady 
gas flows [1^ [27! and analytic solutions for continuum 
scale problems |28l - l3H . However, to the authors knowl- 
edge, for liquids at micro-scales the research has been 
limited to steady flows. To this end, channel flows with 
oscillatory wall movement, the so-called Stokes second 
problem, appears as a natural extension of steady Cou- 
ette flows investigated by Thompson and Troian. Such a 
flow can be encountered iir several microsystem applica- 
tions, such as microaccelerometers, iirertial sensors, and 
resonant filters [3^ . 

In this paper the effects of unsteady flow on slip in sim- 
ple liquids at the solid interface are presented. The nu- 
merical experiments conducted indicate that slip velocity 
is also dependent on fluid acceleration, iir addition to the 
shear rate of the fluid. Previous studies on unsteady flows 



in microchannels by Khalcd and Vafai [g^l , and Matthews 
and Hill [2^ have used the Navier slip model, while 
Bahukudumbi et al. [s^ used the Maxwell slip model for 
their analysis. Both Navier and Maxwell's model sug- 
gest the dependence of slip only on shear rate. To con- 
firm our hypothesis, Maxwell's theory for slip in rarefied 
gases for a steady fiow is extended to an unsteady 
case. Although the original derivation is inspired based 
on some characteristics of gases, we will verify here that 
an analogous formulation is also valid for liquid flows. 
The unsteady slip model developed shows the existence 
of an additional acceleration term. Using the momen- 
tum equation the acceleration term can be rewritten as 
gradient of shear rate. Hence showing that slip is depen- 
dent on shear rate and its gradient. Scaling the model by 
wall velocity and characteristic length for Stokes second 
problem, collapses the data to give a universal boundary 
condition for unsteady and steady flows. It is seen that at 
the limiting case when unsteady flow tends to steady fiow 
the model reduces back to Navier and Maxwell's model. 
We also introduce a non-dimensional number that helps 
in explaining the transition of slip boundary condition 
from finite slip to perfect slip and determining when this 
would arise. Furthermore, the occurrence of hysteresis 
in unsteady flows is shown. Hysteresis is observed when 
comparing slip velocity with shear rate and the fluid ve- 
locity with wall velocity. 

Details of numerical experiments and code validation 
are specified in sections [TT] and IHII In section IIVI un- 
steady slip model is derived. Non-dimensionalization of 
the model is discussed in section [V] Results and conclu- 
sion are made in section |VT] and |VlIl respectively. 



II. NUMERICAL SETUP 

The molecular dynamic simulations presented in this 
paper are performed using the LAMMPS package [s^ ]. 
The problem geometry used is similar to that of Stokes 
second problem which is achieved by selecting the height 
of the fiuid channel to be greater than the Stokes pene- 
tration depth as shown in Fig. [T] The penetration depth 
as determined from the analyti cal s olution of the Stokes 
problem, is given hy 6 = 6.51-\/ j^/w, where v is the kine- 
matic viscosity of the fluid and uj is the frequency of wall 
oscillation. The channel dimensions and number of wall 
and fluid atoms for different wall densities and various 
time periods are presented in Table ID The fluid's ini- 
tial state is modeled as a face centered cubic (fee) struc- 
ture with the X direction of the channel being aligned 
along the [112] orientation of the face-centered cubic lat- 
tice. The wall is comprised of two layers of atoms oriented 
along the (111) plane of fee lattice. The wall atoms are 
flxed to their lattice sites. The bottom wall is kept sta- 
tionary while the top wall is subjected to an oscillatory 
motion, defined by 

X = Asin{ujt), (1) 
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FIG. 1. Schematic of the problem geometry, where h is the 
height of the fluid channel and U is the wall velocity. 
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TABLE 1. Dimensions of the fluid domain and the number of 
fluid and wall atoms are enlisted for different values of time 
period. Time period T, wall number density pm, length x, 
height y, and width z of the fluid channel, number of fluid 
N f and wall A^u, atoms are the variables mentioned. Here r is 
characteristic time of the Lennard- Jones potential and p is the 
fluid number density. The minor changes in the dimensions 
with change in wall density are done to make the problem 
geometry symmetric and periodic. 

where A is the amplitude of wall oscillation. Periodic 
boundary conditions are imposed along the x and z di- 
rections. 

The pairwise interaction of atoms separated by a dis- 
tance r is modeled by the Lennard Jones potential 




where e and cr arc the characteristic energy and length 
scales. The cutoff radius, Tc is 2.2cr where the potential 
is zero for r > r^. 

The fluid is maintained in its equilibrium state hav- 
ing a number density p — O.Slcr"^ and temperature 
T ~ lAks/s- The temperature is regulated by a thermo- 
stat which simulates the transfer of heat from the system 
to an external reservoir. A Langevin thermostat with a 
damping coefficient of F = 1.0r~^, where r = -^/mo^/e, 
is used here. The damping term is only applied to the z 
direction to avoid biasing the flow. The equation of mo- 
tion of the fluid atom of mass m along the z component 
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is therefore given as follows 



mzi 



dVij 



dzi 



mTzi 



(3) 



Here X^j^^i denotes sum over all interactions with i and rji 
is a Gaussian distributed random force. The value of dy- 
namic viscoisty used for the calculations is /i = 0.2eTa~^. 
The equations of motion were integrated using the Verlet 
algorithm [H, [s^ with a time step Tc = 0.002t. 

The simulation is initially run for a time of ^ 600t 
allowing the flow to equilibrate. After which, an ensem- 
ble average of required variables are taken in addition to 
spatial averaging. The spatial averaging is done along 
the length and width of the channel, with a bin height of 
0.25(7. 



III. VALIDATION OF SIMULATION 
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by its limiting value, Lg, versus the shear rate which is 
normalized by the critical shear rate, jc- This primar- 
ily ends up collapsing the data onto one curve given by 
Ls = L°{1 — 7/7c)~^^^. The universal curve shows that 
for a given shear rate the non-dimensionalized slip length 
is independent of the fluid wall properties of the problem 
being considered. 



IV. UNSTEADY SLIP MODEL 

The two widely used model namely Navier's and 
Maxwell's slip model are in essence the same. Both as- 
sume that the slip velocity at the fluid-solid interface is 
proportional to the shear rate of the fluid. Even though 
Maxwell's model js], 113, HI] was established for rarefied 
gases we illustrate here that an analogous formulation 
can be made for slip in liquids. 

Maxwell's theory states that the reflection of fluid 
atoms after colliding with wall atoms can be categorized 
into two types. The two types being specular reflection 
and diffusive reflection as shown in Fig. [31 




Incident fluid atom 



Specularly reflected 



FIG. 3. Schematic of diffusively and specularly reflected fluid 
atom. 



FIG. 2. Universal curve describing the flow boundary con- 
dition. Lb and 7 is scaled by and jc respectively. The 
curve flt is given by Ls = L%{1 — 'y/'jc)"^^'^ in agreement by 
Thompson and Troian's results [^. 

Before going ahead with our present experiment we val- 
idated our code by reproducing Thompson and Troian's 
results . The problem consisted of a Couette flow ge- 
ometry. Simulation was run for different cases having 
varying wall-fluid properties. For each case, slip length 
was computed for increasing shear rate of fluid which in 
turn is governed by the wall velocity. The slip length 
is calculated from the linear Navier boundary condition. 
Ay = Lsj, which can be simplified for Couette fiow to 
{U/j — h)/2 0, where AV is the slip velocity, Lg is the 
slip length and 7 is the fluid shear rate. The shear rate 
is computed as the slope of the velocity profile. The re- 
produced results of Thompson and Troian Q is shown in 
Fig. [2] Considerable agreement was found with their re- 
sults. The non linear dependence of slip length with shear 
rate is illustrated. Also we were able to duplicate the uni- 
versal curve which is a plot of slip length normalized 



Diffusive reflection (u£^y. In this case the incident 
fluid atoms exhibit no-slip and the velocity of the atom 
attained after collision is the same as that of the velocity 
of the wall during collision. This can be imagined as fluid 
atoms being adsorbed by the wall and then put back into 
the fluid in a random direction. Hence, the aggregate 
velocity of the atom after collision being equal to the 
velocity of wall atom at the instant of collision (relative 
velocity with wall being zero) 

UlHtt) = U-{Q. (4) 

Here tc is the instantaneous time of collision and JJ™ is 
the wall velocity. 

Specular reflection (U^"^) : The atoms undergo perfect 
slip and the velocity after collision is the same as that 
before collision 

Ul^itt) = uLit-), (5) 

where Ul, is the fluid velocity at the wall. Further- 
more, the velocity of the atom before collision with wall 
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is deemed to have been obtained by collision with fluid 
atoms located at a distance equal to A from the wall. 
This being the location of the first fluid layer and is 
the time it takes for the fluid atom to travel this distance. 
Here we make an assumption that during the motion of 
the fluid atom in reference it is subjected to a net zero 
force during its journey from the first fluid layer to the 
wall. Thereby, it has a constant velocity 



(6) 



Performing a Taylor series expansion in space at the wall 
the velocity at a distance A from it can be written as 



Ulitc 



,) = Ul (U, - T,) - A 



dy 



(7) 



Discrctizing in time and ignoring all second and higher 
order terms we obtain 



ui{t-)^ui{t,-T,) = ui{t,y 



dt 



A- 



dy 



(8) 



Now, the average velocity after collision Ul,{t'^) is 
given by, 

Ulitt) = aM\tt) + (1 - Ul\tt), (9) 

where ad is called the tangential momentum accommo- 
dation coefficient (TMAC), which gives the fraction of 
atoms undergoing diffusive reflection. Substituting from 
equations ([4]) and ([5]) in the above equation we obtain 

Ul{tt) = cjaU^ {Q + (1 - <Jd) Ulit-). (10) 

The net instantaneous fluid velocity at the wall is given 
as the mean of the velocity before and after collision 



Ulitc) 



ul{t+) + ul{t-) 



(11) 



By substituting Ul{t+) from dTU]), one obtains 

2Ul{tc) = adU^ {Q + (2 - <Td) Ul(t-). (12) 
Now U£ {t~ ) could be replaced from equation ^ to give 
2Ul, (tc) = aaU^ (U) + 2Ul (t,) - + 



(2 - aa) 



dUl , dUf (tc) 



dt 



dy 



(13) 



Rearranging this equation results in 



itc)-Ul (ie) 



(2 - ad) 



dUlitc) , ,dUf{tc) 



dt 



+ A- 



dy 



Slip velocity is given as. 



(14) 
(15) 



Using the definition of slip velocity equation (|14l) can be 
written as 



Us{tc 



(2 ~ ad) 



dUlit,.) , ,dUf{tc) 



dt 



+ A- 



dy 



(16) 



Finally writing dU /dy in terms of shear rate, 7 we get 



Us {tc) 



(2 - ad) 



Cd 



dUl (t,) 
' dt 



A7 



(17) 



Hence, we see that slip velocity has an additional de- 
pendence on the fluid acceleration at the wall in case of 
unsteady flows. 

The simplified Navier Stokes equation is used in order 
to rewrite the acceleration term in the above expression 
in terms of shear rate. Writing down the 2-D Momentum 
equation in x, with no external force or pressure gradient 



du 
'dt 



d 



dx dy 



dx dy 



In addition since the velocity of fluid in the x direction is 
uniform, the rate of change of velocity along x direction 
goes to zero. Hence the Navier Stokes equation reduces 
to, 

- ■ (19) 



du 
'dt 



dy" 



Now considering that, 7 
equation as, 



dy 



one can rewrite the above 



dul 

dt 



.9 fdUl\ 
dy\dy) 



97 
dy 



(20) 



Substituting this in the equation for slip velocity, equa- 
tion ([T7)) . results in 



Us (tc) 



(2 - ad) 
ad 



(21) 



For steady flow, slip velocity is seen to be proportional 
to shear rate of the fluid as hypothesized by Navier and 
Maxwell. But here, we show that for an unsteady prob- 
lem we have an additional contribution from the gradient 
of shear rate. Hence it is seen that slip velocity for an 
unsteady case is proportional to the linear sum of shear 
rate and gradient of shear rate of the fluid at the wall. 

In the limitng case when the flow approaches steady 
state and the acceleration goes to zero the unsteady slip 
model reduces back to that given by Maxwell, 

'2-(t/ 



Us = 



A7 



(22) 



Comparing this to Navier's slip bondary condition 

Us = Ls 7L . (23) 



one can write the slip length, Ls ~ 



X. A ad value 



of one can be compared to a no-slip boundary condition 
in liquids where there is no relative velocity between the 
wall and fluid. While ad equal to zero corresponds to a 
boundary condition exhibiting perfect slip. 
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V. NON-DIMENSIONALIZING THE 
UNSTEADY SLIP MODEL 

Non-dimensionalization of the slip model is performed 
with an attempt to simplify and parametrize the equa- 
tion. It is also aimed at identifying a non-dimensional 
parameter that could characterize unsteady slip. The 
unsteady slip model established in previous section is 
scaled by the wall velocity, Uw, and the length scale 
associated with Stokes second problem, ^Ivjuj. Non- 
dimcsnioalizing using this velocity and length scales, one 
obtains 



2_ad 



A dm 




(24) 

where (.)* represents non-dimensional quantities. Notice 
that, Tc = ^, where c is the speed of sound in the fluid 
at a given temperature. Therefore 



u: 



o-d 



Grid 



d^U*f dU} 
dy*"^ dy* 



(25) 



This could be also written as. 



u: 



l: 



d^U] dU 



wher 



(^) 



Qy*2 Qy* 



A* and Cnd = 



(26) 



/2uui 



, . . „ „„ -T^^ is the non- 

o-d J 2c 

dimensional parameter that controls the unsteady term 
based on the type of flow. One can see that for a limiting 
case of w = 0, which corresponds to a steady flow prob- 
lem, the equation reduces back to Maxwell's slip equa- 
tion. 



VI. RESULTS 

The primary information that is extracted from the 
simulations is the fluid velocity. The following steps are 
taken before the data from the simulations is used for any 
analysis. First, a reference plane is defined at a distance 
of 0.5(7™-^ away from the wall lattice site. This is the loca- 
tion at which the fluid variables at the wall is computed. 
Secondly in order to obtain a well resolved velocity pro- 
file the Levenberg-Marquardt method [H, H^] is used to 
fit the analytical solution of the Stokes second problem 
to the data. This is done so as to limit the noise resulting 
from taking the derivatives of the velocity profile. Differ- 
ent cases of wall-fiuid properties considered are enlisted 
in Table HH 

Thompson and Troian Q observed that slip length for 
fiuid shear rates much lower than the critical shear rate 
is constant. But for shear rates in the vicinity of the 
critical value it becomes highly nonlinear As we are 
here considering an unsteady problem the various fluid 
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TABLE II. Four different cases with varying wall-fluid prop- 



erties were considered, e™^ and a™^ are the Lennard- Jones 
parameters for fluid-wall interaction. 



and wall variables are time varying, and so is slip length 
if the shear rate is in the vicinity of its critical value. 
Hence, the slip length corresponding to each run is calcu- 
lated by fitting the data using the Lavenberg-Marquardt 
method [H, Ho] and the slip model given in equation [26l 



A. Verification of slip model 

Several numerical simulations were conducted in or- 
der to verify the unsteady slip model established in the 
previous section. The variation of slip length with time 
period, shown in Fig. SJ^a), gives us considerable insight 
into the unsteady problem. It is seen that the slip length 
approaches infinity as the wall oscillation frequency tends 
to infinity. This implies that the wall is oscillating at such 
a high frequency that no information can be passed on 
by the wall atoms to the underlying fluid atoms, thereby 
giving a perfect slip. As we increase the time period we 
observe that the slip length decreases and approaches an 
asymptotic finite value which corresponds to those shown 
by Thompson and Troian in the steady state Couette flow 
problem [9[. 

In Fig. m^b) the non-dimensionalized slip length, L* 
also mentioned as is plotted against x which is given 

as Cnd ^y'i + where a; is a sum of the parameters 

that influence unsteady slip. The data for different cases 
of wall-fluid properties and different wall velocities col- 
lapse onto a single curve. This curve is well described by 
the form y = a/x^ where the dashed line represents the 
values a = 0.5 and h = 1.0. The exact fitting coefficients 
to our experimental data were calculated to be a = 0.59 
and h = 0.96, given by solid line which closely matches 
the dashed line. Hence, it is seen here that similar to 
Thompson and Troian proper scaling leads to a universal 
curve which holds for unsteady and steady flows. With 
the help of this curve one can find the slip length for a 
given value of shear rate and its gradient. 



B. Relevance of Non-dimensional number 

In the non-dimensional number Cnd, \phAjj is the phase 
velocity of the propagating velocity profile in y-direction 
whereas c is the speed of sound through the medium. The 
ratio can be physically seen as the ratio of propagation 
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(a) 



C. Effects of acceleration on slip 
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FIG. 4. (a) Slip length versus time period. As the unsteady 
problem tends to a steady problem, steady state value of slip 
length is recovered, (b) Non-dimensional slip length versus 
sum of the non-dimensional shear rate and the product of 
non-dimensional number C„d and gradient of shear rate. Here 

y = L*g and x = Cnd '^gyJi + ■ The curve fitting parameters 
are a = 0.59 and h — 0.96. It is seen that all the data for 
various cases considered collapse onto a single curve given by 
the fit. 



of momentum through the medium to the speed of prop- 
agation of sound through the medium. The maximum 
speed at which momentum can be transferred through 
a medium is at the speed of sound. Hence, when the 
phase speed is of the order of the speed of sound for the 
medium, the fluid will exhibit perfect slip as the wall 
would not be able to transfer momentum on to the fluid. 
The number can be used to determine when the bound- 
ary condition would change to perfect slip from a finite 
or no-slip condition. 



In order to analyze the affect of acceleration on slip 
a wall-fluid property is chosen which is shown to ex- 
hibit no-slip at the fluid-solid interface for steady flow 
by Thompson and Troian This corresponds to case 
1 in Table [TTl Wall oscillation having an amplitude of 
IOct and time period of 40t is chosen such that the shear 
rate is considerably below its critical value. In Fig. [SI 
comparison of the velocity profile obtained from simula- 
tion is made with Stokes analytical solution with no slip 
boundary condition. This done for an instantaneous time 
t = 40t for which distinct slip at the wall is observed. 




-Analytic with no slip 
■MD simulation 



-0.5 



0.5 



u/U 



1.5 



FIG. 5. Velocity profile obtained from simulation with Stokes 
analytical solution for no-slip boundary condition is com- 
pared. The profile is for case 1 of wall-fiuid properties and 
for an instant of time equal to the time period of wall oscilla- 
tion. A distinct presence of slip at the wall is seen. 



1. Phase lag of fluid velocity due to wall acceleration 

The presence of fluid slip at the wall indicates the ex- 
istence of a lag of fluid velocity with respect to the wall. 
This is investigated by comparing fluid velocity for vary- 
ing time periods of wall oscillation with the wall velocity 
as shown in Fig. [51 This is done to see the effects of ac- 
celeration on the fluid velocity. The comparison is made 
for case 1 of wall-fluid properties. The amplitude of wall 
oscillation is calculated such that the amplitude of wall 
velocity remains constant. It is seen that the phase lag 
reduces and the amplitude of fluid velocity increases with 
increasing time period relative to the wall. An increase in 
time period results in a decrease in fluid inertia. Hence, 
it can be said that as we approach a steady state, by 
increasing time period, we recover the no-slip boundary 
condition for case 1 . Leading to the conclusion that iner- 
tia affects the slip of fluid at the wall. Thereby, conflrm- 
ing the general slip model derived in previous section. 
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time/time period 

FIG. 6. Comparison of wall velocity and fluid velocities at 
time periods 40r and 200t for a cycle of wall oscillation having 
wall- fluid properties corresponding to case 1. An amplitude 
of 10a and 50a corresponding to the time periods of 40r and 
200t are chosen in order to achieve the same amplitude of 
wall velocity in both cases. An increase in wall-fluid phase 
lag and a decrease in the amplitude of fluid adjacent to the 
wall are observed by decreasing time period. 

2. Phase lag of fluid velocity due to wall hydrophobicity 

The variation of fluid velocity with time for different 
cases mentioned in Table |TT] is shown in Fig. [71 This is 
done for a fixed time period of 40t and fluid acceleration 
of 10a. The purpose being to see how the fluid velocity 
behaves with varying hydrophobicity. The fluid velocities 
lag with respect to the wall velocity. There is an increase 
in phase lag with increasing hydrophobicity. Also, a re- 
duction in the amplitude of the fluid velocity is noticed. 
The phase lag and amplitude reduction are a result of 
the increase in slip. Tang et al [2^ also observed similar 
phase lag of fluid velocity for different Stokes numbers 
and TMAC in their Lattice Boltzmann simulation of os- 
cillatory gas flows. 



D. An explanation for the occurrence of slip 

Increase in hydrophobicity leads to an increase in slip 
at the wall, which has been observed by several re- 
searchers for the steady case @-[l3|- Strictly speaking, 
the no-slip boundary condition is only valid if the flow 
adjacent to a solid surface is in thermodynamic equilib- 
rium [2^, [3l|. For fluid flow in small-scale systems, the 
collision frequency is not high enough to ensure thermo- 
dynamic equilibrium, thus a certain degree of tangential 
velocity slip must be allowed [2^. Also, these collisions 
should occur during a time interval smaller than that of 
the smallest time scale for flow changes. Harris and Rice 
calculated the relaxation time for liquid argon at 90 K as 
0.5r while for gaseous argon at 300 K to be of the order 
of lOOr [4l| . The relaxation time required for bulk liquid 
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time (r) 

FIG. 7. Velocities of wall and adjacent fluid in one wall os- 
cillation cycle. The wall oscillation has an amplitude of lOcr 
and time period 40r. An increase in wall- fluid phase lag and 
a decrease in the amplitude of fluid adjacent to the wall are 
observed by increasing the wall-fluid hydrophobicity. 

argon is small in comparison to the time scale of molecu- 
lar collisions. However, as shown in Fig. [5] fluid layering 
is observed close to the wall. This leads to the reduc- 
tion of liquid density to the order of gas density adjacent 
to the wall. This reduces the number of fluid atoms in- 
teracting and undergoing momentum transfer with the 
wall thereby increasing the required relaxation time [2^ . 
Hence, the fluid atoms adjacent to the wall do not have 
sufficient time to equilibrate and the transfer of momen- 
tum from the wall is only partial; therefore, resulting in 
slip. 

55| 1 1 1 




0.5 1 1.5 

density(l/cr^) 



FIG. 8. Variation of fluid density along the height of fluid 
channel. Layering of fluid is observed close to the wall. Red 
dashed line marks the reference plane O.Scr'"-'^ away from the 
wall lattice site. The plot corresponds to case 1 with an am- 
plitude of lOcr and 40r. As the layering is symmetric about 
X — z plane only the top half is shown. The ordering is preva- 
lent up to a distance of 5 — 6cr from the reference plane with 
decreasing amplitude beyond which the bulk density of the 
fluid is obtained. 
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E. Boundary condition for an unsteady flow 



The main aim of the studies pertaining to slip length 
is to find the exact boundary condition, which is slip 
velocity at the wall. We find that similar to slip length, 
slip velocity also collapses onto a single universal curve 
when non-dimensioanlized with the appropriate scaling 
parameters as shown in Fig. HI The scaled slip velocity 



(a) 




FIG. 9. Data for different cases collapses onto a single curve 



au'f 



defined hy y = ^, where y ^ Us and x = C„d-gpi- + 
The parameters for the curve fit are a=f .05 and b=f .f4. 



is defined as, 



Us = u: 



(27) 



where U* is the slip velocity scaled by velocity scale of Uw 
as mentioned in section [V] The plot data uses the max- 
imum slip velocity and the instantaneous value of fluid 
shear rate and gradient of shear rate at the wall corre- 
sponding to the occurrence of maximum slip. A univer- 
sal curve is well described by the form y = a/x^ where 

y = Us and x = Cnd + f-nd the dashed line rep- 
resents the value a = 1.0 and h = 1.0. The exact fitting 
coefficients to our experimental data were calculated to 
be a = 1.05 and b = 1.14. 

Proper scaling leads to the collapse of data which sug- 
gests a universal boundary condition for any wall-fiuid 
property. The universal curve equips us with a very po- 
tent tool in determining the slip at the solid-fluid bound- 
ary for a unsteady and steady flows having the knowledge 
of fluid shear rate. Thereby eliminating the need to per- 
form extensive molecular dynamic simulations which are 
computationally expensive. 
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FIG. 10. Slip velocity as a function of (a) wall velocity and (b) 
fluid shear rate exhibiting hysteresis is shown. The plots are 
for case 1 and an amplitude of lOcr and 50(7 corresponding 
to a time period of 40r and 200r was chosen in order to 
maintain the same maximum wall velocity. All the curves are 
for a complete cycle of oscillations. It is seen that the area 
confined by the hysteresis loop decreases with increase in time 
period. 



F. Hysteresis 

The lagging of fluid velocity due to fluid inertia sug- 
gests the presence of hysteresis and forms the motivation 
to explore it in unsteady oscillating flows. In Fig. [TOla) 
the slip velocity at the wall is plotted against the wall ve- 
locity and in Fig. [TUfb) change in slip velocity is plotted 
against the fluid shear rate. Two different time period 
of wall oscillation, 40t and 200r, are considered in each 
figure to study the effects of inertia on hysteresis. A hys- 
teresis loop is formed in both the figures and the area 
confined within the loop is seen to decrease with increas- 
ing time period. 

The effect of inertia in developing hysteresis as seen in 
Fig. [TO] is explained here. This could be described by vi- 
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sualizing fluid layers in the vicinity of the wall during the 
time when the wall is approaching an extreme location 
in the cycle and then changing direction. As the wall ap- 
proaches the extreme point in the oscillation all the fluid 
atoms are moving in the same direction as the wall. But, 
when the wall changes its direction the fluid adjacent to 
the wall does not change direction instantaneously as a 
result of inertia of fluid atoms in this and the neighboring 
layers. This leads to a change in the magnitude of fluid 
velocity, depending on the direction of motion of wall, 
although they are evaluated at the same location of wall 
in the cycle. 

While there is a distinct loop observed in Fig. [TOla') at 
a time period of 200r the loop in Fig.fTOlb) collapses into 
a single curve. The area formed by slip velocity and wall 
velocity has the dimensions of energy per unit mass while 
that formed by slip velocity and fluid shear rate has the 
dimensions of acceleration. Hence, the area in Fig. [TUTa) 
could be possibly interpreted as the additional amount of 
energy the wall needs to transfer to the fluid below it in 
order to attain thermal equilibrium. This can be related 
to the lack of thermal equilibrium of fluid atoms adjacent 
to the wall which leads to slip. The area in Fig. [TUfb) 
having the dimensions of acceleration can be correlated 
to inertia in the fluid. As the time period is increased 
from 40r to 200t we essentially reduce the inertia in the 
flow thereby getting closer to steady state which leads to 
the collapse of the loop. 

The hysteresis here can be seen to be analogous to 
the commonly known magnetic hysteresis [4^ . When 
a ferromagnetic material is subjected to magnetic field 
the atomic dipoles aligns itself with the field but when 
the field is reduced to zero, partial alignment is retained. 
In order to demagnetize it a field must be applied in 
opposite direction. This results in the difference in the 
path taken by the two legs in the hysteresis loop. 

VII. CONCLUSION 

A series of numerical experiments with different wall- 
fiuid interaction properties and varying amplitudes and 
time periods were performed for a time-periodic oscilla- 
tory Couette fiow problem. An increase in slip is ob- 
served as compared to the steady Couette flow problem. 
This increased slip was attributed to fluid inertial forces 
not represented in a steady fiow. A case having wall- 



fiuid properties that showed no-slip for Thompson and 
Troian's steady Couette flow experiment is chosen. For 
this case when the fluid is subjected to an oscillatory flow 
a distinct slip is observed which confirms the increase in 
slip in an unsteady fiow. To gain a deeper insight into 
the cause of the increased slip an unsteady slip model is 
established based on the Maxwell's slip boundary con- 
dition. The dependence of slip on acceleration in ad- 
dition to shear rate is shown. By writing acceleration 
in terms of shear rate, it is shown that slip at the wall 
depends on gradient of shear rate and the shear rate of 
fiuid at the wall. For the limiting case of steady fiow the 
model reduces back to the Maxwell's model. This pro- 
vides a more accurate prediction of slip for unsteady flow 
problem rather than simply using the steady Navier's or 
Maxwell's slip model. Non-dimensionalizing the model 
by scaling the problem by the wall velocity and the char- 
acteristic length of Stokes second problem leads to the 
collapse of data onto a single universal curve. Thereby 
with the knowledge of shear rate of fluid, slip length can 
be calculated using the universal curve without having to 
perform any computationally expensive MD simulations. 
A key non-dimensional number, defined as the ratio of 
phase speed to speed of sound, helps in explaining and 
characterizing the transition of slip boundary condition 
from finite to perfect slip is also identified from this non- 
dimensionalization. Phase lag in fluid velocity relative to 
wall is observed. The lag increases with decreasing time 
period of wall oscillation and increasing hydrophobicity. 

Hysteresis is observed while comparing slip velocity 
with wall velocity and shear rate. The cause for hys- 
teresis can be attributed to the inertia of fiuid. It is seen 
that the area formed by the loop decreases with increase 
in time period of wall oscillation. For the case of slip ve- 
locity versus the shear rate the loop collapses to a curve 
as the flow tends to a steady flow at an increased time 
period. In the case of slip velocity versus the wall veloc- 
ity the area inside the hysteresis loop can be related to 
the loss of energy transfer from the wall to the fiuid. For 
the loop formed by slip velocity and shear rate the area 
can be said to be equivalent to inertia of the flow. 
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